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One of the striking features of particle production at high beam energies is the near equal abun- 
dance of matter and antimatter in the central rapidity region. In this paper we study how this 
symmetry is reached as the beam energy is increased. In particular, we quantify explicitly the en- 
ergy dependence of the approach to matter/antimatter symmetry in proton-proton and in heavy-ion 
collisions. Expectations are presented also for the production of more complex forms of antimatter 
like antihypernuclei. 

PACS numbers: 25.75.-q, 25.75.Dw, 13.85.Ni 

Keywords: Production of antimatter, relativistic heavy-ion collisions 



Ph. 

CD 



> 
ON 

rn 
o 



X 



I. INTRODUCTION 

One of the striking features of particle production at 
high energies is the nearly equal abundance of matter 
and antimatter in the central rapidity region P, 0- It 
is believed that a similar symmetry existed in the ini- 
tial stage of the universe. It remains a mystery how this 
symmetry got lost in the evolution of the universe reach- 
ing a stage with no visible amounts of antimatter being 
present. Closely related to the matter-antimatter sym- 
metry is the production of light antinuclei, hypernuclei 
and antihypernuclei at high energies, especially in view 
of the recent observation of the anti ^He nucleus by the 
STAR Collaboration Q. 

Since the first observation of hypernuclei in 1952 Q 
there has been a steady interest in searching for new 
hypernuclei and exploring the hyperon-nucleon interac- 
tion which is relevant (see e.g. 0, Ig) for nuclear physics. 
Hypernuclei decay with lifetime which depends on the 
strength of the hyperon-nucleon interaction. While sev- 
eral hypernuclei have been found since the first observa- 
tion no antihypernucleus has ever been observed until the 
recent discovery of the antihypertriton in Au-t-Au colli- 
sions at v^sTra = 200 GcV by the STAR Collaboration 
at RHIC [3]- The yield of (anti)hypernuclei measured 
by STAR is very large, in particular they seem to be 
produced with a similar yield as other (anti)nuclei like 
e.g. (anti)^IIe. This abundance is much higher than mea- 
sured for hypernuclei and nuclei at lower energies Q . It is 
of interest to understand the nature of this enhancement, 
and for this, the systematics of antimatter production 
in high energy hadron-hadron and heavy-ion collisions 
should be investigated. 

The analysis of particle production assessing the de- 
gree of thermalization of the particle source has been un- 



dertaken since many decades [ol-fT^. It has been found 
that the thermalization assumption applies successfully 
to hadrons produced in a large number of particle and nu- 
clear reactions at different energies |15l - [l7| . This fact al- 
lowed to estimate thermal parameters characterizing the 
particle source for each colliding system which is relevant 
for the understanding of the thermal properties of dense 
and hot matter and for studies of QCD phase transitions 

El [a. 

In this paper, using the energy dependence of thermal 
parameters obtained from the statistical thermal model 
analysis of particle yields in heavy-ion collisions [lO, HH 
we present the model estimates for the (anti)hypernuclei 
multiplicity that can be directly compared to the recent 
results obtained in central Au-|-Au collisions at RHIC. 
We discuss systematics of (anti)matter production at dif- 
ferent energies. We also make predictions of (anti)matter 
and (anti) hypernuclei production at the Large Hadron 
Collider (LHC). 

Recently, a very interesting analysis of the production 
of light nuclei, hypernuclei and their antiparticles in cen- 
tral heavy-ion collisions was performed in Ref. [2^ within 
the statistical thermal model. It was found that ratios 
of hypernuclei to Lambda exhibit an energy dependence 
similar to the K^/tt^ ratio with a clear maximum at low 
energy. The present work is considered to be an extension 
of the analysis performed in Ref. ^2^] . 

Firstly, we compare the statistical thermal model re- 
sults on (anti)baryon production in heavy-ion and in 
proton-proton collisions. This, in general, requires the 
knowledge of the energy dependence of thermal param- 
eters in p-fp collisions which are proposed in this pa- 
per based on the recent data. In this context, we study 
quantitatively how the matter/antimatter symmetry is 
reached as the beam energy is increased both for p-f-p and 
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heavy-ion collisions. We also discuss the role of the 
strangeness content of particles and quantify different an- 
tibaryon/baryon ratios in p+p and in heavy- ion collisions 
at SPS, RHIC and LHC energies. 

Secondly, we compare predictions of the statistical 
thermal and coalescence models for different ratios of 
(anti)nuclei and (anti)hypernuclei in the context of re- 
cent STAR data obtained in central Au-t-Au collisions at 
the top RHIC energy. 

The paper is organized as follows: In section II we 
discuss features of the statistical thermal model. In 
section III we compare the antibaryon/baryon ratios 
in p-l-p and heavy-ion collisions and obtain the energy 
dependence of thermal parameters in p+p collisions. 
We demonstrate the scaling behaviour of the an- 
tibaryon/baryon ratio with their strangeness content. 
In section IV we apply the thermal and the coalescence 
models to the production of nuclei and hypernuclei 
and their antiparticles. We also make predictions for 
(anti)nuclei and (anti)hypernuclei yield ratios at LHC 
energy. In section V we summarize our results. 



II. THE STATISTICAL THERMAL MODEL 

The statistical thermal model assumes that in a high 
energy collision at freeze-out all hadrons follow equilib- 
rium distributions. The conditions at chemical freeze-out 
where inelastic collisions cease are given by the hadron 
abundances, while the particle spectra offer insight into 
the conditions at thermal freeze-out where elastic col- 
lisions cease. Once thermal parameters are fixed, the 
hadron gas partition function gives all primordial ther- 
modynamic observables of the system. The exact form 
of the partition function, however, depends on the statis- 
tical ensemble under consideration. 

Within the grand-canonical (GC) ensemble, the quan- 
tum numbers of the system arc conserved on average 
through the action of chemical potentials [l^. In other 
words, the baryon _B, strangeness S and the charge con- 
tent Q are fixed on average by the fin, fJ-s and fig chem- 
ical potentials respectively. For each chemical potential 
one can introduce the corresponding fugacity A = e^^'^ 
where T is the temperature of the system. 

In the GC ensemble the density of hadron species i 
with the mass mi, the quantum numbers Bi, Si and Qi 
and with the spin-isospin degeneracy factor gi is express 
through the second order modified Bessel function K2{x) 
as 

n,(T,/.B,/.5,/iQ) = ^mM'AfA§'X2(^). (1) 

The above form, valid in the Boltzmann approximation, 
is easily generalized to the quantum statistics [1^ [2l| . 

In the application of the statistical thermal model the 
chemical potentials /is and /xq are typically constrained 
in the initial stage by the strangeness neutrality condi- 
tion and by the fixed baryon-to-charge ratio. However, it 



is well established, that the usual form of the statistical 
thermal model formulated in the GC ensemble cannot 
be used when either the temperature or the volume pa- 
rameter V or both are small [l^, [2^ . As a thumb rule 
one needs VT^ > 1 for a grand canonical description 
to hold [23 - [26| . This condition is not usually justified 
in p+p collisions, requiring canonical (C) formulation of 
strangeness conservation. The exact strangeness conser- 
vation causes a suppression in particle ratios of strange 
(or multi-strange) hadrons to pions or any strangeness 
neutral particles as compared to the corresponding ratio 
in the grand canonical limit. The key parameter gov- 
erning this effect can be quantified by the strangeness 
correlation volume UM ■ 



HI. PRODUCTION OF ANTIBARYONS 

The chemical freeze-out conditions in heavy-ion colli- 
sions at various energies were shown to follow a curve 
in the temperature-baryo-chemical potential plane (27j 
which has been phenomenologically parameterized as 

T{fiB)^a-b^i%-cn% (2) 

Hb = d/{l + c^snn) (3) 

with a = 0.166 ± 0.002 GeV, b = 0.139 ± 0.016 GeV^^ 
c = 0.053 ± 0.021 GeV-3, d = 1.308 ± 0.028 GeV and 
e = 0.273 ± 0.008 GeV~^. This parametrization is quan- 
titatively similar to the one proposed in Ref. [l^ and 
results in a very satisfactory description of different parti- 
cle excitation functions measured in nucleus-nucleus col- 
lisions. 

The increase of the antimatter to matter ratio in 
heavy-ion collisions with the center-of-mass energy has 
been observed by the NA49 [13, [Ml and the STAR ^ 
Collaboration. Fig. [1] shows changes of the p/p ratio 
with collision energy at mid-rapidity in central heavy-ion 
and in p+p collisions. The data from NA49 and STAR 
Collaboration are compared with new results from the 
ALICE Collaboration There is a clear increase of 
this ratio towards unity, indicating approximate symme- 
try of matter and antimatter at the LHC energy. There 
is also a clear increase of the p/p ratio when going from 
heavy-ion towards p+p collisions. 

In Fig. [1] data are compared with statistical thermal 
model results. In heavy-ion collisions these model calcu- 
lations were done using the energy dependence of model 
parameters as described by Eqs. ^ and ([3]). There is 
a clear agreement of model predictions with data. For 
p+p collisions no systematic analysis of model parame- 
ters with energy were performed till now. 

If (anti)nucleons are directly originating from a ther- 
mal source, then from Eq. ([T]) (i.e. neglecting feed-down 
from resonances) it is obvious, that the p/p densities ra- 
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exph2/XB/T], 



(4) 



is entirely quantified by the ^b/T value. Thus, an in- 
crease in the p/p ratio from heavy-ion to p+p collisions, 
seen in Fig. [U is due to a decrease in the ^mb/T value. 

To extract the corresponding hb and T at fixed energy 
in p-l-p collisions we have used the THERMUS code [20, 
[2]| which correctly accounts for feeding corrections to 
(anti)nucleons from decays of heavier resonances. 

The p/p ratios measured in p-fp collisions, shown in 
Fig. [TJ have been fitted using the statistical thermal 
model by varying only the parameters d and e in Eqs. ^ 
and ([3]). We have used the same T{^s^^) dependence 
for p-l-p as for heavy-ion collisions. This is justified by 
the observation that at high energies there is no notice- 
able change in T between central and peripheral heavy- 
ion as well as p-l-p collisions [33|. The resulting baryo- 



chemical potential hb is shown in the lower part of Fig. [T] 
by filled circles. In addition, applying the parametriza- 
tion of ^b(^snn) as in Eq. ^ we have found that the 
parameters corresponding to p+p collisions are 

(5) 



MS = rfpp/ (1 + Cpp^/snn) 

with dpp = 0.4 GeV and epp = 0.1599 GeV-i. The 
solid line in the lower part of Fig. [T] represents the energy 
dependence of /is in p-l-p collisions obtained with the 
above parameters. For comparison also shown in this 
figure is the energy dependence of the value of /.ts in 
heavy-ion collisions. From Fig. [T] it is clear that at mid- 
rapidity, the ^B is always lower in p+p than in heavy-ion 
collisions. This observation reflect the fact that at mid- 
rapidity the stopping power in p+p collisions is less than 
in heavy-ion reactions. The change of p/p ratio with 
energy in p+p collisions is quantified in the upper part 
of Fig. [l]using parametrization of /is (-y/sNN) adjusted for 
p+p collisions. 

For baryons carrying a number of Ns (anti)strange 
quarks the antibaryon/baryon ratio (again neglecting for 
the moment feed-down from resonances): 



TIB 



exp[-2(AiB-iV5Ai5)/T], 



(6) 



is modified by the strange chemical pot ential. As /xs is 
always smaller than ^b (see e.g. [TgI. |34||). the above ra- 
tios should appear ordered with the strangeness quantum 
numbers, i.e. the higher iVg, the smaller the difference 
between antibaryon and baryon. 

Figure [2] shows data on strange antibaryon/baryon ra- 
tios from the SPS and RHIC energies and their compar- 
isons with the model calculations using the THERMUS 
code. The data and the model results both in p+p and 
heavy- ion collisions are in good agreement. There are 
clear trends in strange antibaryon/baryon ratios already 
expected from the simplified Eq. ([6]). (i) With increasing 
strangeness quantum number, the antibaryon/baryon ra- 
tios are increasing and approaching to unity, (ii) Heavy- 
ion collisions exhibit smaller B/B ratios as compared to 
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FIG. 1: The p/p ratio (upper part) and the corresponding 
/i_B (lower part) as a function of -^snn are shown. The solid 
circles are results from p+p collisions and the open squares are 
results from heavy-ion collisions [ll.B [30l - [3^ . The dashed line 
is the parametrization for heavy-ion collisions from Ref. [2^ 
while the solid line is the new parametrization for p+p colli- 



p+p collisions due to different [Ib values as mentioned 
before. This is well seen at SPS energies, where the 
difference in [Lb in p+p and Pb+Pb systems is larger 
than at RHIC. (iii) The differences between heavy- ion 
and p+p collisions decrease with increasing y^snn. At 
LHC energies the p/p ratio is close to unity and there- 
fore, the abundances of strange baryons are roughly as 
large as those of their antiparticles. 



IV. PRODUCTION OF (ANTI)NUCLEI AND 
(ANTI)HYPERNUCLEI 

A. Comparison to data from RHIC 

The production of antimatter compared to matter, 
being expressed by ratios of antibaryon/baryon yields, 
was shown in the last section to be well described by 
the statistical thermal model. Thus, it is of interest to 
verify whether the recently observed production of light 
(anti)nuclei including (anti)hypertritons (^H) in heavy- 
ion collisions at RHIC by STAR Collaboration Q also fol- 
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FIG. 2: Antibaryon to baryon ratios sorted according to their 
strangeness content. Circles (solid horizontal line) refer to 
p+p collisions data (model calculations) and open squares 
(dashed horizontal line) refer to heavy-ion collisions data 
(model calculations). The upper part shows results at the 
SPS and the lower part at the RHIC energy. 



lows a pattern expected in the statistical thermal model. 

Studying the antinuclei/nuclei ratio in the statistical 
thermal model an extra factor of /xb is picked up each 
time the baryon number is increased. Thus, each nuclcon 
adds a factor of fiB in the exponent of the Boltzmann 
factor in Eq. ([T]) . The production of nuclear fragments is 
therefore very sensitive to the value of the baryo-chemical 
potential and thus could be used for a precise determi- 
nation oi ^B- 

The deuterium has an additional neutron and the an- 
tideuterium/deuterium ratio in the statistical thermal 
model is given by 



^ =exp[-4Mi3/r], 
rid 



(7) 



which then is '-^(p/p)^. 

If the nuclei carry strangeness, this leads to an ex- 
tra term and the ratio of antihypertriton/hypertriton 
reads 



713 u 



(9) 



In mixed ratios, i.e. using ratios of different nuclei (or 
antinuclei), there appears an extra factor due to differ- 
ent degeneracy and masses, e.g. in the statistical thermal 
model 



9lU {m^^nf ^2(m3H/r) 

53Hc ("l3Hc)^ iC2(m3Hc/r) 



exphA^s/r]. (10) 



Figure [3] shows comparisons of the statistical thermal 
model results on different (anti)nuclei ratios with recent 
experime ntal data from the STAR Collaboration. These 
■^He and ■^He yields have been corrected for contamina- 
tion from hypertriton and antihypertriton decays assum- 
ing a decay branch ratio of 25% and consequently in the 
model such decays have not been included. 

In the statistical thermal model, following Eqs. ([5]) and 
(|9|). ratios of (anti)nuclci/nuclei are entirely quantified 
by the Hb/T and Hs/T values. From Fig. [3] it clear that 
using the thermal parameters at chemical freeze-out ob- 
tained from the analysis of particle yields at RH IC, t here 
is an excellent description of measured ratios of ■^He/'^He 
and ^H/^H. However, deviations are seen on the level of 

mixed ratios, ^H/^He and ^H/^He. 

In elementary collisions nuclei and antinuclei as well 
as hypernuclei and antihypernuclei can be produced by 
direct pair production. In heavy-ion collisions, due to fi- 
nal state correlations, a different production mechanism 
opens up through hadron coalescence. Indeed, produc- 
tion of nuclei in Pb+Pb collisions at ■y/SNN= 17.3 GeV at 
CERN SPS [s^ have been found to be consistent with a 
coalescence picture, while this was not the case in p+Be 
collisions at the same energy. 

In the most straightforward coalescence picture the ra- 
tios of different (anti)nuclei can be directly related to 
ratios of hadronic yields. In particular. 



^He ^ ppn ^ p^3 
^He ppn p 



3h 
All 

A^^ 



pnA 
pnA 



(11) 



(12) 



thus should be similar to the square of the antipro- 
ton/proton ratio if decay contributions of heavier res- 
onances to nucleon yields are neglected. The ^He has 
three nucleons and the corresponding anti-^He/^He ratio 
is given by 



713 Ho 



exp[-6/XB/r], 



(8) 



and 
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(14) 
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FIG. 3: Comparison of results from the STAR Collaboration 
with the statistical thermal and the coalescence model. For 
the latter both experimental values [s^jls^l (dashed lines) and 
values from the statistical thermal model (dash-dotted lines) 
have been used. 



Comparing results of the statistical thermal model 
with the coalescence framework, one sees that some ra- 
tios are quite similar. Indeed, from Eqs. ^ and (|TT|) as 
well as from Eqs. ([9]) and ([T2|) it is clear that neglecting 
feed-down from resonance decays the statistical thermal 
model coincides with coalescence predictions on the level 
of 3lfe/3He and pI/iH ratios (see also_Ref. [H). Thus, 
as long as the key input ratios p/p and A/ A are in agree- 
ment with a thermal descriptions, the measured ratios do 
not allow to distinguish the two mechanisms. However, 
differences betw een these models are seen on the level 
of mixed ratios, ^H/^He and \ll/^Iie, due to different 
masses of nuclei. From Eqs. ((T0| and ([TS]) one finds that 
when neglecting binding energy of nuclei and feed-down 
corrections the statistical thermal model differs from the 
coalescence framework by a factor of (l/3-|-2mp/3mA)^^^. 
Consequently, the statistical thermal model results for 
^H/^He and ^H/'^He ratios (solid lines in Fig. [3]) are 
lower than those obtained in the coalescence picture us- 
ing the (anti)A/p ratios from THERMUS. 

The results from the coalescence model [H, [s^ are 
compared to data from the STAR Collaboration and the 
statistical thermal model predictions in Fig. [3l The coa- 
lescence estimate has been done using the p/p, A/A, A/p 
and A/pratios both measured by the STAR Collabora- 
tion 0,l36(, Jil] (dashed lines) and from the THERMUS 
calculations (dash-dotted lines). 

We note that in coalescence picture the equilibrium 
abundances of particle yields are not required. Conse- 
quently, (anti)nuclei produced from the off-equilibrium 
medium can lead to particle ratios being in agreement 
with the simple coalescence estimate discussed above. 
However, this is not anymore the case for statistical ther- 
mal model which requires statistical order of particle 
yields in the final state. 



B. Predictions for RHIC and LHC 

In the previous section we concentrated on the statisti- 
cal thermal model description of (anti)matter production 
in heavy-ion collisions up to RHIC energies. In the fol- 
lowing we extend our discussion to higher incident ener- 
gies and quantify differences between p+p and heavy-ion 
collisions. 

In Fig. 21 we compare p+p and heavy- ion collisions at 
■y/sNN = 200 GeV. In the context of the statistical ther- 
mal model the difference between these two colliding sys- 
tems is caused by different values of /is and by the effect 
of canonical suppression in p-f-p collisions. The ratios of 
antinuclei-to-nuclei without strangeness content are only 
affected by the baryo-chemical potential which at mid- 
rapidity is smaller in p+p than in heavy-ion collisions as 
discussed earlier. 

Figure HI upper part, nicely demonstrates that with 
increasing mass the effect of fiB becomes stronger, yet, a 
strangeness content causes an opposite trend as discussed 
earlier. The ratio of hypertriton-to-"^He and the corre- 
sponding antimatter ratio show the effect of the canoni- 
cal suppression reducing the yield of (anti)baryons carry- 
ing strangeness. For the chosen correlation volume with 
Rc = 1.5 fm the difference is not dramatic but very no- 
ticeable. 

The effect of increasing collision energy is demon- 
strated in the middle part of Fig. 2] Here, the differ- 
ences between the antimatter/matter ratios in heavy- 
ion collisions at RHIC and LHC are essentially due 
to the decreasing value of /i^. At LHC, the chemi- 
cal potential is smaller than 1 MeV resulting in the 
antimatter/matter ratio being close to unity. The ra- 
tios of the (anti)hypernuclei/'^(anti)He remain nearly un- 
changed from RHIC to LHC since here the effect of hb 
is only due to the strange chemical potential which is 
small. These ratios are dominated by mass differences 
and degeneracy factors. 

The expectations for LHC energies are studied in more 
detail in the lower part of Fig. 21 by comparing the 
p+p and Pb+Pb collisions. For simplicity, in both cases 
the collision energy of 7 TeV has been chosen. The ratios 
do not change between ^snn = 2.76 TeV and 7 TeV. The 
antimatter/matter ratios are hardly changed from p+p to 
heavy- ion collisions. All antiparticle/particle ratios are 
close to unity. The ratios of (anti)hypernuclei/^(anti)He 
exhibit the influence of the canonical suppression for 
the correlation volume (see Section |lT]) corresponding to 
Rc = 1.5 fm [2^. For larger Rc the canonical effect is 
reduced and already for i?c = 4 fm is hardly visible. 

The predictions of the statistical thermal model for ra- 
tios of particles with different masses are shown in Fig. [S] 
The calculations have been performed for Au+Au colli- 
sions at y'snn = 200 GeV. 

For the results presented in the preceding figures, we 
used the freeze-out temperature according to Eqs. ^ 
and ([3]). It is clear that ratios of nuclei with different 
masses are strongly influenced by the value of the freeze- 
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FIG. 4: Comparison of different particle ratios calculated in 
the statistical thermal model using T — 170 MeV. Upper fig- 
ure: for p+p and heavy-ion collisions at y^snn — 200 GeV. 
Middle figure: For heavy-ion collisions at different collision 
energies. Lower figure: Prediction for p+p and Pb+Pb colli- 
sions at x/snn = 7 TeV. 



out temperature. Figure O displays statistical thermal 
model results obtained with T varying between 110 MeV 
and 170 MeV. These calculations are compared with the 
recently measured value from the STAR Collaboration 
including the observation of anti-alpha particles (si. l4Ci|. 
More data are needed before a freeze-out temperature for 



antinuclei can be concluded. 
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with the recently measured values [J HO]. 



V. SUMMARY 

We have discussed in a systematic manner the produc- 
tion of (anti) matter in high energy collisions within the 
statistical thermal model. We have focused on a general 
comparison of the production of antibaryons and antin- 
uclei. The variation of the p/p ratios with ^snn be- 
ing different for p+p and heavy-ion collisions has been 
used to obtain the parametrization of the energy depen- 
dence of thermal parameters in p+p collisions beyond the 
SPS energy. We have demonstrated the scaling behavior 
of the (anti)baryon/baryon ratios with the strangeness 
quantum number and the changes in these ratios between 
p+p and heavy-ion collisions with y/suN- 

We have compared the measured ratios of nuclear and 
anti-nuclear fragments in heavy-ion collisions with the 
statistical thermal model and with the coalescence con- 
cept. Based on the successful description of existing data, 
we have made predictions for (anti)mattcr production in 
p+p and heavy-ion collisions at LHC energies. 
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